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Observational evidence suggests that our universe is presently dominated by a dark en- 
ergy component and undergoing accelerated expansion. We recently introduced a model, 
motivated by string theory for short-distance physics, for explaining dark energy without 
appealing to any fine-tuning. The idea of the transplanckian dark energy (TDE) was based 
on the freeze-out mechanism of the ultralow frequency modes, u(k) of very short distances, 
by the expansion of the background universe, uj(k) < H. In this paper we address the issue 
of the stress-energy tensor for the nonlinear short- distance physics and explain the need to 
modify Einstein equations in this regime. From the modified Einstein equations we then 
derive the equation of state for the TDE model, which has the distinctive feature of being 
continually time-dependent. The explanation of the coincidence puzzle relies entirely on 
the intrinsic time-evolution of the TDE equation of state. 
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1 Introduction 



Cosmological observations of large scale structure, SNla, age of the universe and cosmic 
microwave background (CMB) data, strongly indicate that the universe is dominated by a 
dark energy component with negative pressure |I[ . Besides the difficulty of coming up with 
a natural explanation for the smallness of the observed dark energy, an equal challenge is 
the "cosmic coincidence" problem. 

Recently we proposed a model [§] for explaining the observed dark energy without 
appealing to fine-tuning or anthropic arguments. This model is based on the nonlinear 
behavior of transplanckian metric perturbation modes which was motivated by closed string 
theory |]3], |[] and quantum gravity ||. The transplanckian dark energy (TDE) model was 
based on the freeze-out mechanism of the short- distance modes with ultralow energy, by 
the expansion of the background universe, H, and it naturally explained the smallness of 
the observed dark energy. 

In this paper we study the stress-energy tensor of the TDE model in order to calculate 
the equation of state for these short- distance stringy modes. As we will show, the frozen 
tail modes start having a negative pressure of the same order as their positive energy 
density soon after the matter domination era. Thus it is only at low redshifts that they 
become important for driving the universe into an accelerated expansion and dominate the 
Hubble expansion rate H . A distinctive feature of the TDE model is that its equation of 
state, wh, is always strongly time-dependent at any epoch in the evolution of the universe, 
(e.g. wh = —1/3 during radiation dominated era but it becomes wh = —1/2 at matter 
domination). It becomes more and more negative at late times until it approaches the 
limiting value wh = — 1, after the matter domination time, t eq . 

The calculation of the components of the stress-energy tensor, T^ u , namely the pressure 
and energy density, is given in Section 2. Due to the nonlinearity of short distance physics, 
Bianchi identity is generically violated for all these models. Therefore one needs to modify 
the Einstein equations, (T^), such that the modified ones satisfy Bianchi identity. 

From physical considerations, the need for modifying Einstein equations in the nonlinear 
regime of short distance physics is to be expected, due to nonequilibrum dynamics of 
the short distance modes. In practical terms this is not an easy fair to carry out in an 
unambiguous way, for a simple reason: we do not have a unique effective theory valid at 
transplanckian energies or a lagrangian description of the theory in this regime |J. The 
only information available to most transplanckian models [0, |S|, [It], [II], [12] is the field 
equation of motion (with a few exceptions like Jl3|). Nevertheless all these models do violate 
Bianchi identity and the energy conservation law, if T^ u is not modified accordingly. 

Based on the equation of motion as our sole information for short-distance physics, we 
therefore use a kinetic theory approach for modifying Einstein equations in the absence 
of an effective lagrangian description. The assumption made is that a kinetic theory de- 
scription of the cosmological fluid is valid even in the transplanckian regime. Despite its 
nonlinear behavior at short distances, this imperfect fluid shares the same symmetries, 
namely homogeneity and isotropy, as the background Friedman-Robertson- Walker (FRW) 
universe. Then the corrections r^ u to the stress energy tensor T^ v will also be of a diagonal 
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form JTJ] 



r^u = (e + U)u^u u + Ylg^ u 



(1) 



In Section 3 we explore the observational consequences of the model with the puzzle 
of 'cosmic coincidence' in mind. A summary is given in Section 4. A discussion of the 
nonequilibrum dynamics and distribution function for the transplanckain modes, as well as 
details of averaging of their energy and pressure, are attached in the Appendix. 



The Equation of State from the Modified Einstein 
Equations 



2.1 Analytical expression for T, 



Transplanckian models that investigate the sensitivity of the CMB spectrum or Hawking 
radiation to short-distance physics, all introduce a nonlinear, time-dependent frequency for 
the very short wavelength modes 0, |8|, 



<"\p]=f\p] = f[k/a] 



(2) 



The physical momentum p is related to the comoving wavenumber k by p = k/a, with a the 
scale factor. Most of these models lack a lagrangian description and, all the information 
they propose about short- distance physics is contained in the mode equation of motionf} 



[n+Lu(k,a] 



0. 



(3) 



The expectation that Einstein equations will not hold unless they are modified in the 
nonlinear regime of short distance physics is fully reasonable and it is based on the fact 
that Bianchi identity and energy conservation law will be violated due to the nonlinear time 
dependence of u. In terms of kinetic theory, the time-dependence of the group velocity v g 
indicates departure from equilibrium [H| (see Appendix). Here we study the modifications 



of T^v for a specific class, the TDE model || . Our approach is based on kinetic theory and 
the pressure modifications are obtained through balance equations. 

In the TDE model we are considering, the dispersed frequency for short- distance metric 
perturbation modes is: 



o; 2 [p] = p 2 £[p/p c 



a 



V 

1 + u 
exp[2p/p c 



+ 



e 3 u 



u 



(4) 
(5) 
(6) 



lr The up- term collectively denotes the generalized frequency that appears as a mass squared term in 
the equation. Depending in the particular problem studied it may also include other terms to it, like 
for example, the coupling of the modes to the curvature of the universe, a" /a, if the equation under 
consideration is that of metric perturbations. 
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where p c is of order the Planck mass or string scale M, p is the physical momentum, e« 
arbitrary constants. The maximum of u[p] is around p ~ p c . The frequency function 
behaves as 

u 2 [p] « p 2 (l + 0(p 2 /M 2 )) p<M (7) 

for the modes in the sub-planckian regime, and like 

uJ 2 [p\ ~ \/ei + e 3 p 2 exp[-2p/M] p>M. (8) 

for those modes in the TP regime. The nonlinear exact function Eq (|5| for the frequency 
can be fitted to a good accuracy to oj[p} 2 ~ C osh[p/p -i] 2 • 

Let us refer to the wavepackets of the modes centered around a momentum pi as par- 
ticles. Then, their group velocity v g = du/dp is time dependent through its nonlinear 
p-dependence, and is different form the phase velocity, v c = uj/p. 

Lorentz invariance is broken due to the nonlinearity at short distances. Therefore, the 
fixed cuttoff scale p c = M, together with all the transplanckian modes pick a preferred frame, 
the CMB frame. This frame is freely falling along the comoving geodesies, with respect 
to the physical FRW Universe?]. Sometimes we will refer to transplanckian modes as the 
modes inside a small box with fixed Planck size, l p = 1/M, in the preferred frame, since their 
wavelength Xpp < l p is smaller that the "size of the box" , p > M. In this picture, Lorentz 
invariance is broken in the small box but restored in the large box with size L = a/M, i.e 
the Universe. Thus the physical momenta modes for the "small box" bound observers in 
the preferred frame are the comoving wavenumber modes for the "outside" observers, in 
the Lorentz invariant FRW Universe, that " see" the preferred frame in a free fall. 

Let us address below the issue of how the energy density components behave with 
time, prior to the pressure modifications. The short distance modes are out of thermal 
equilibrium, due to their nonlinear frequency and group velocity, v g ^ 1. Meanwhile a 
thermal state is restored at large scales, (A ^> lp), where the frequency is nearly linear and 
thus v g ~ 1. Thus we need to average the contribution of the short distance modes to the 
energy and pressure in the Universe, over many of their wavelengths, in order to obtain an 
effective large scale thermal state. That is why in obtaining the equation of state, < Wi > 
prior to the pressure modifications Ilj, for the transplanckian modes, the averaging is done 
in time-scales of cosmological order. Details of averaging are provided in the Appendix, 
5.2. The equation of state < Wi > prior to viscous pressure modifications, is obtained from 
the expression < Wi >= with < pi >= — < pi > — < (a/3)dpi/da > and a the scale 
factor. Based on the behavior of v g with p, we divide the dispersion function into 4 regions 
(see Fig.l) 

• Region : Linear regime, up to p ~ p c , such that uo[p\ ~ p. These modes behave as 
radiation, Eq.(^), with the averaged pressure expression being p = po/3. 

• Region I : Around the maximum of the dispersion function, up to some value pp > M 
in physical momentum, where u can be expanded in a polynomial series, and the leading 

2 See for a very nice treatment of issues related to a fixed physical cuttoff in a preferred frame. 
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order terms in Eq. (|]) are the first ones. Region I is dual to region 0. 

oo[p] « p (oq + ai (p/M) + a 2 (p/M) 2 + ■) , (9) 
where a« are constants, (a 2 < 0). We use Eq. @ to estimate the energy density 

" = l£ , ' <i " 3 ( a ° + ai ss + 02 (^ (10) 

CM 4 fa 0t 4 a a 5 



a* VT (l «- 1) + ^- 1) + '"J (U) 

M 4 

oc — , (12) 
a 4 

where = k B /M > 1, = 0(1). The constant C = \j3 p \ 2 denotes the Bogolubov 

coefficient squared, which in our model does not depend on the momentum p || . Therefore, 
pi behaves like radiation plus 0(1/ a 2 ) corrections in its averaged equation of state, (pi) = 
(1/3 — A/a 2 + • Regions I and contribute to the radiation energy component in 

the Universe. 
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• Region II: From some mode, ps 3> M, onwards denned such that its frequencies can 

be best fitted to an exponential dependence on p, u[p] ~ pexp[—p/M]. The energy for this 
region is 

p n ~ — J k dkk 3 exp[-2k/aM] (13) 

= C^-(F[x B ]-F[x H ]) , (14) 

where, 

F[xi] = d + M + ^1 + 3 ) exp[-2 Xl /a] . (15) 



and Xi = ki/M. Since xb > 1 then ~ ^x%exp[— 2xb/o]- Thus p/j behaves as matter 

when averaged over many oscillations and its averaged pressureQis (p n ) = (—B/a)(pn) ~ 0. 
Also since xh > %b > 1 then F^b] >> F[xh\- 

• Region "H" : This is our "tail"H, defined as the part of the graph for which the 
frequency of the modes is smaller than the Hubble parameter, H. The functional behavior 
of the frequency with p is the same as in region II, therefore the averaged pressure expression 
for this region is the same as that of Region //, that is (pn) — 0. But the lower limit of 
integration kn (or pu) is given by the physical condition of the freeze-out of the modes by 
the expansion of the background universe, 

u H \p H ]=H. (16) 

This region includes the modes from pu to oo in the range where u is exponentially sup- 
pressed. The energy density of the tail is 

pH ~ 7^ exp[ ~ 2kH/aM] (17) 

Notice that due to the freeze-out, the evolution of the fc# mode is highly nontrivial and 
thus corrections to the averaged pressure term (pn) — will be important. 

Modes in the tail, between pn to oo, behave differently from the other modes, since 
their time-dependence is controlled by the Hubble expansion, Eq. (|16"D. On the other 
hand, all modes with momenta p < Ph redshift in the same way with the scale factor, 
towards decreasing values, i.e the linear regime^. Nevertheless, these regions, (0, /, II), also 
receive small modifications to their pressure term from the deep transplanckian regime. We 
show below that the modifications due to the p#-defrosting effect, are non negligible and 
important only in the highly nonlinear regime, around pn- 

Now, we would like to estimate the corrections to pressure, n^, for all these regions, 
with the notation, P 4 for the effective modified pressure 



Pi -> + n 



3 See Appendix for details of averaging. 

4 Modes in the linear regime are referred to as " normal modes" . 
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where the index runs to % = 0, /, II, H. The averaged unmodified pressure expressions, (p%), 
are 



(p ,i)^(l-A/a 2 + ---)(p 0tI ), (19) 
(j) II>H )~{-B/a + ---){p II>H ). (20) 

The terms that go as inverse powers of a can be neglected and A,B are numerical con- 
stants related to the averaging (see Appendix for details). We refer to these expressions as 
"bare pressures" in order to distinguish them from the viscous pressure modifications terms 
(defined below). 



In a similar manner to particle creation cases |L6[ in imperfect fluids (17], jig] , the highly 
nontrivial time-dependence of the mode pn and the transfer of energy between regions, due 
to the defrosting of this mode across the boundary p H , gives rise to pressure corrections in 
the fluid energy conservation law. The defrosting of the modes results in a time-dependent 
"particle number" for regions near p^- From kinetic theory we know that this "particle 
creation", (the defrosting of the modes), gives rise to effective viscous pressure modifications 
TA [13. The term IT denotes the effective viscous pressure modification to the "bare" 



pressure, (pi). 

The criteria we will use for modifying is that Bianchi identity must be satisfied [|19j 
with the new expressions for pressure^], Pi, 

+ 3H{ Pi +Pi + IL,)] = Zi[ Pl + 3H{ Pi + P^] = , (21) 

with i = 0,1, II, H. Let us write this expression explicitly in terms of its energy and 
bare pressure components, and collect the contributions of regions 0,1 into one combined 
radiation energy, p R = po + pr- 

pii + 3H{p II +p II )+p R + 3H{p R + p R ) = -3HU n , (22) 

p H + 3H{p H +p H ) = -3HTl H , (23) 

where pn,H — 0, p R — 1/3. So, we have imperfect fluids in regions II and "H" , and 
eventually their energy is transferred, due to the redshifting effect, to regions and I, which 
is why these regions also receive pressure modifications. Nevertheless, the viscous pressure 
corrections to the "radiation" modes are very small since the energy and the volume of 
phase space occupied by them is very large (a 3 times larger the Planck size volume). These 
regions are in a nearly equilibrium situations, (see Appendix). All pressure corrections are 
estimated below. 

Let us find out Ilj, in order to solve Eqs. ( ]22| ) and (^). As explained, the presence 
of Hi is due to the exchange of energy between the two regions, from the defrosting of the 
modes pu at the boundary. This is directly related to the time dependence of the boundary 
Ph, which in turn is going to be controlled by the Hubble parameter H. In essence, there 
is an exchange of modes between region (R + //) and "H" . Although the specific number 



5 From here on we drop the (• • ■) notation and denote the averaged "bare" pressure by simply pi instead 
of (ft). 
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of particles^] in each of these regions, Nu and Nh, is not conserved, their rate of change, 
in the physical FRW Universe, is related through the conservation of the total number of 
particles which contains both of these components 

N T = , (24) 

Each component satisfies[] 



Nu = rjjNn, (25) 
N H = (3H-T H )N H , (26) 

where the "decay rates" of the regions r$ account for the rate of change in the number of 
their "particles" (modes), due to the defrosting effect. 

The system is not yet in equilibrium. The change in the number of "particles" is giving 
rise to the effective viscous pressure, IL> Even prior to the freeze-out effects, that is, even 
for Tjjji = 0, the short distance modes in region II and III were out of thermal equilibrium, 
due to their nonlinear frequency and group velocity, v g ^ 1. 

The contribution terms to pressure, IT, are related to Tj through [jlT|, 



3HU H = -{p H +p H )T H , (27) 
3HU H = -[(pii+pii) + (pR + pR))rn. (28) 

Therefore, Eq. (|23| ) reduces to 

p H +{3H -T H ){p H +p H ) = 0, (29) 

which can be also recast as: 

Ph = —(ph+Ph), (30) 
n H 

where nu = Nh / a 3 is the "particle" number density for the region of modes from pn to 
infinity. The flow of particles is described by nn = nnu a , with u a the unit 4-velocity vector 
of the fluid. Notice that since the group velocity in the H-region is is negative, particles in 
this region flow in a direction v g which opposite to the direction of their momenta, k. The 
rate r#, calculated below, is positive, and the increase in the number of particles Nh as 
given by Th does not allow the energy density of the "tail" to redshift as fast as matter. 
This indicates that although small, pu eventually will come to dominate the total energy 
density. 

6 We are loosely using the term particle here to refer to the wavepackets of the transplanckian modes, 
centered around a momenta pi 

7 Vector objects related to the flow direction of the fluid are denoted in bold letters, e.g. N; = NiU a with 
u a being the unit 4-velocity of the fluid and the corresponding modulus of this vector being iVj. Notice 
that the factor (3HNh) in Eq. (|2^) is related to the fact that the preferred frame for the tail modes is 
falling along comoving geodesies in the Universe. 



7 



The number of "particles" Nh contained in the tail regime, in its preferred frame, is 
given by 

/•CO 

N H ~ C H dkk 2 exp[-k/aM] 

J kg 

~ C H (aM)k 2 H exp[-k H /aM}. (31) 

where Ch is the constant proportional to the Bogolubov coefficient^ (3 p . We can now 
calculate the energy transfer, due to the defrosting of the modes kn, between the tail 
region and region // from the balance equation for Nh, Eq. (|26|), where Nh is: 



r°° f k \ 

N H ^ Ch J k dkk 2 i—j exp[-k/aM] - C H k H exp[-k H /aM]kH 
~ 3HN H -C H k 2 H exp[-k H /aM](k H - Hk H ) 

In the last line we have used the approximation in Eq. (ftTD, and: 

Ph Ph/M .H. 



p H I-Ph/M'H^ 
derived from Eq. flTB|). When pn 3> M (which always holds), we have for Tr'- 



(33) 



TH^m+^-^m 1 T otal )> ( 34 ) 

si Z 

where w to tai = Ptotall Ptotal- Therefore, when w tota i -> -1 then IV reaches its limit, T H -> 
3iT. can not change anymore once this limit is reached because the Hubble constant 
and the mode pn freeze to a time-independent value. Notice that Th is positive for all 
equations of state Wtotai < 1 and thus it slows down the dilution of the tail with the scale 
factor. 

We can repeat the same procedure for the modes in region 0, J, // in order to obtain a 
closed equation for prji , similar to Eq. (p9f), i.e. that is given entirely in terms of pnji 
and p RJI 

Pii,r + (3if - T n ) {p H ,R + Pii,r) = , (35) 

where we have used Eq. (p8|) . 

Let us now try to relate Tr to Tjj. The total number of not frozen particles, Njj, in 
the region from zero to kn is given by 

N n ~ C H [4:M 3 + Mphujr] ■ (36) 



8 In our model |/3 P | 2 calculated in M resulted in a scale invariant spectrum. It does not depend on 
the wavenumber p, the reason why it can be pulled out of the integral and factored into the coefficient 
C'h = N/3p, with N an overall normalization constant we are keeping for the sake of generality. 
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From the total balance equation for the particle number between the two regions, (R + 
II) and region "H" in the comoving volume, we have Ntotai = 0, where Nn = IV/Nii and 
N H = Mp 2 H Exp(-p H /M)u = Mp H Hu. Thus 

r " = iJf^H ^- 3H) - (37) 

and N T = N n - C H (Mp H u H ) = AC H M 3 . In obt aining the scalar quantity for the number 
of particles Nt from their flow Nt , the negative sign picked up in the second term in Nt 
is related to the fact that the flow of the tail's defrosted modes is in the opposite direction 
to their momenta, due to their negative group velocity. Therefore, by plugging in the 
expression of Nu from Eq. (^), we get that in the limit Nu ^> MpuH, Fu is smaller 
than Fh and negative, given by the expression: 

F n = -{3H-F H )^-^-. (38) 

Since PhH 2 << M 2 H then y = '^j^- ~ 0(H/M) is going to be much less than 1 for as long 
as the expansion is not dominated by the tail. From the condition ujh[ph] = H, and the 
time evolution of the physical momentum pu in Eq. (|33|), we have that pujpu —* —H/H, 
when ph 3> M. The exact value of y does not matter and it is small. The tail domination 
case, when pn becomes comparable to pu, should be treated separately since Tjj — > 0. The 
pressure modification II// increases the dilution of pu as determined by the equation for 
Pu. This equation shows that due to the modified pressure effects, pu goes to zero faster 
than a matter energy density component. 



2.2 Equation of State 

In this part we calculate the effective equation of state for all the regions, from the pressure 
expressions, pi, Hi, that were obtained in the previous section. Starting with region "H", 
we have 

= (T H -3H)(l + w H ) = -^(l + w total )(l + w H ) (39) 
Ph * 

where we have defined 

H 3 

- = --(l + Wt otal), (40) 

with Wtotai referring to the effective equation of state for the total energy density. The 
"effective" equation of state for the tail can be read from this expression to be 

l + w H = -(1 + w tota i) (l + w H ), (41) 
with wh — 0. The time evolution for pu is 

3 r 

Ph = pn(0)exp[-- J (1 + w tota i)(l + w H )d\na\ (42) 
= p H (0)exp[-3 f (1 + w H )d\na] . (43) 



9 



During radiation domination, Wtotai = 1/3, then the effective equation of state for the tail 
is Wh = —1/3; for matter domination, w to tai = 0, then % = —1/2; at the start of the 
accelerated expansion, q = 0, pn = p#, we have w to tai = —1/3, % = —2/3; and finally, if 
the tail dominates, then the only solution to the Friedman equation is given by w tota i — wh, 
with wh = —1. Therefore, the tail starts behaving as dark energy only recently, when its 
equation of state wh becomes close to the limiting value, % = — 1. 

By the same procedure, we can now estimate the effective equation of state for pn in 
terms of V u : 

= T n -3H = -3H[1 + ( 1 ~ Wtotal )y] (44) 
Pn 2 

Thus the effective equation of state for region II, obtained from Eq. is: 

l + w„ = (l + y l ' W 2 total ) (1 + wn) , (45) 

where y = ) = 0(H/M) and w n ~ 0. 

The time evolution of the pn energy density component is 

Pn = pn(p)exp[—3 J{1 + w n )d\na] , (46) 

In a radiation dominated universe, H = M/a 2 therefore (H/M) = (1/a 2 ), and matter 
dominated, H = M/a 3 ^ 2 , so then (H/M) = (1/a 3 / 2 ). The point is that the correction y 
to the matter equation of state for region 77 is really small for up to the equality time. 
During most of the history of the Universe, y goes as an inverse power of the scale factor 
a. So region II behaves pretty much like matter. The special era when the tail eventually 
dominates the expansion and H becomes a constant, (at pn = constant, wh = w to tai — — 1) 
is discussed below in Sect. 3. 

Similarly, by repeating the same steps, from Eq. (|22"D , it can be shown that the modifi- 
cations to pressure for pn, regions 0,1, are very small indeed. Thus their effective equation 
of state remains very nearly that of radiation, wr = 1/3. In order to avoid repetition, we 
will not carry out the calculation for the effective equation of state of pn, as the procedure 
is essentially the same as for pji, and it results in an effective radiation equation of state 

l+w R =(l + y l ~ W 2 total ){l + w R ) , (47) 



3 The Issue of Coincidence and Comparison to Ob- 
servation 

From the computation of the pressures pn and IT, it is clear that the initial radiation is 
redshifted faster than the other components of the total energy density. We can ask the 
question at what time, t eq , the pn components of matter becomes comparable to radiation 

PR*Pii**f- = \KM 2 , (48) 
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with 

pu ~ Cp 3 B exp[-2p B /M} ~ ^ e 2 g M 2 . (49) 
From the above equation we obtain at t eq : 



M 

Vb ^ y In 



2C ^ (50) 



3HLM 2 



On the other hand, from Eq. (jig) we have: 

p# ~ Mln 



(51) 



if 

and comparing Eqs. (0) and flST]) , it is clear than PB(t eq ) < Puiteq), and therefore 

PH^eq) < Pll(teq) ■ (52) 

So matter-radiation equality takes place well before the eventual "tail" domination. From 
the equations of state, wu and wu, Eqs. ( filD and (f45|), we have that pjj always dilutes 
faster than pjj. Thus the inequality in Eq. (60) holds true not only at t — t eq but at 
all earlier times before t eq . Generally, there may be other sources of matter and radiation 
in the Universe, besides the contribution from the transplanckian modes. Although these 
components would not be affected by the viscous pressure corrections, IL, their contribution 
should be including in the Friedman equation when determining the equality time, a(t eq ) = 
a eq . Since their effect to the expansion is well studied and known, here we chose to focus 
our attention only on the role of the transplanckian modes. 

Let us now estimate the time at which the tail takes over to dominate the expansion 
and address the issue of the cosmic coincidence. As we will see below, the effective equation 
of state wu for p u changes from w u ~ to wu > 0. Let us start by asking at what time, 
cide, we have 

Ph = Pu = ptotai/2 (53) 

or in terms of the density parameters Qh = ■ From the Friedmann Eq. for the expansion 
and the relation of wh to w tota i it is straightforward to find out that at a = cl^e we have: 

2 1 

w H {.a DE ) = -- , w to tai(a DE ) = -- , (54) 

and therefore 

, (Q) x 2/3(«w-l) / (0 )\l/2 / o (0)\ 1 /2 

with pf^ being the value of the i-th component at equality time, a eq . It is interesting 
to notice that Wtotai = — 1/3 corresponds to the transition time where the deceleration 
parameter, 

q ~ -(3w tota i + 1) , (56) 



11 



changes sign and goes through zero. This means that acceleration starts at the same time 



as the dominance of the tail, clde, i- e., 



cide- Using the Friedman expansion law, 



we can find the solution for the scale factor a, after the time a q : 



(0) (o) 
Pii , Ph 



H\a q ) 



3M 2 3M 2 



a q y a q 



— + — 

a J a 



(57) 



Therefore, 



a/a q 



a/ a q 



i \ (a/a,) 2 + 1 



d(a/a q 



H{a q 
V2 



<t-t q ) 



(5f 



This integral can be done exactly and it is messy. The important point about it is that it 
gives a power law accelerated expansion, — ~ t n with n > 2. Clearly the tail is behaving 
as dark energy and it is dominating the expansion soon after a,. 

Let us understand physically what is happening around the time a = a q , and why 
a q = clde- As showed in Sect. 2, due to the strong coupling of the tail evolution to the 
Hubble constant, H, and therefore a coupling to ptotai, wh becomes more and more negative 
with decreasing values of w to tai, soon after a eq . Thus the tail starts behaving as dark energy, 
dominates the expansion, and approaches its limiting value wh = — 1 only at late times, 
when other energy contributions to ptotai become negligible. From the Eqs. ([34]) and (|38|), 
when pu ~ pu we obtain Tn ~ — 2yTn- When the tail comes to dominate, the only 
solution to the Friedman equation is T H — > 3 if, and then w to tai = wh — — 1. This means 
that the tail dominates the expansion (w to tai — — 1) very fast and, around that time pu 
has become nearly zero due to the Tn viscous pressure corrections. Recall, that in this 
estimation we assumed that around the time a q , pu is the only source of matter. The fast 
dilution of pn as compared to pn, due to the viscous pressure effects of Tn, is the reason 
why q = occurs at the same time as the tail dominance, a^E- With no other sources of 
matter, w to tai almost immediately goes from w to tai = _ 1/3 to w to tai — — 1- Thus a DE has 
occurred very recently indeed. If we consider other matter contributions in the Friedman 
equation, that are independent of pn, Tj, then the time the expansion takes between the 
start of the accelerated expansion a q , and the time of tail dominance time over ptotai, (i-e. 
when Wtotai = Wh — ~ 1) becomes a bit longer. This time interval from a q to present is 
also the time interval when the tail has acquired a dark energy equation of state, until it 
reaches its limiting value, wh = —I- 

Therefore, cosmic coincidence is explained naturally by the intrinsic time evolution of 
the effective equation of state for the tail, wh ■ 



4 Summary 

Models of nonlinear short distance physics discussed recently in literature ||, |9|, |K], |TT 



12], |13| , |15| , usually introduce a time dependent frequency, at the level of the equation of 
motion for the field. As a result Bianchi identity is generically violated which indicates 
that Einstein equations need to be modified in the high energy regime. It is difficult to 
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do so without an effective lagrangian description of the theory in the transplanckian (TP) 
regime. 

We therefore used a kinetic theory approach, in order to estimate the short distance 
modification to the cosmic fluid stress energy tensor for the model of 0. It is not clear 
to us whether this procedure determines the modifications in a unique unambiguous way, 
or whether fluid idealization and the assumption that kinetic theory remains valid at such 
high wavenumber modes, is a good approximation. Nevertheless, we believe that without 
an effective lagrangian, kinetic theory is the only available tool to get some sensible results 
for the contribution that TP modes have in the long wavelength regime. 

In a previous paper f2|, we showed that the energy contribution of the tail modes is 
comparable in magnitude to the observed dark energy in the universe. In this work we 
calculated the effective equation of state, wh, for these tail modes in order to address the 
cosmic coincidence issue and showed that the tail modes behave as dark energy only at late 
times. 

The tail has an exponentially suppressed frequency and all the modes with uj < H are 
frozen out by the expansion of the background universe. However, due to the short distance 
pressure modification, the tail does not always behave as dark energy. The highly nontrivial 
time-dependence of tail's dominant modep#, tracks the evolution of the total energy density 
ptotai through its strong dependence on the Hubble constant, H. The dependence of pu on 
H is given by the freeze-out condition. 

As a result of the coupling of the pu mode to ptotai, the tail equation of state wh follows 
the evolution of w to tai, such that Wh — ^ Wtot ^ l ~ 1 ' . For this reason, wjj acquires more and 
more negative values as w to tai decreases, from radiation to matter. The tail has a slower 
dilution with the scale factor compared to the other components and it starts dominating the 
expansion and behaving as dark energy only recently, from the time when the deceleration 
parameter q changes sign at time a q . From this point, a q , ( with w to tai = —1/3, wh = —2/3) 
and onwards, the tail drives the Universe into an accelerated expansion, and soon reaches 
its limiting value of w to tai — wh — — 1 with ptotai — Ph- Therefore, cosmic coincidence in 
this model is explained naturally from the time evolution of the tail's % = f[w to tai]- This 
is the most important result of this paper. 

The TDE model was motivated by closed string theory of Brandenberger- Vafa model |3J] . 
Therefore its observational implications may be explored as indirect string signatures. Some 
of the distinctive features of the TDE model are the predictions thatQ: the change in sign of 
the deceleration parameter, q = 0, occurs at the same time as the start of tail dominance, 
ie. the time when the tail energy is half of the total, a q = clde] this accelerated expansion 
occurs very recently; due to the viscous pressure effects, the matter contribution from the 
TP modes goes through a change of its equation of state such that it starts decaying faster 
than normal matter. The latter effect shortens the time the Universe takes to change the 
parameters from w to tai - -1/3, q = to the time when w to tai - -I, Ptotai — Ph- To get 
these numbers, we ignored other matter sources and also the short distance corrections to 
the equations of state (wi) that go as inverse powers of the scale factor, 0{l/a n ). Perhaps 
these corrections may be important recently, in delaying the time it takes the tail to behave 

9 These predictions are in the absence of other matter sources. 
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as dark energy, wh — — 1- These features can be scrutinized with respect to observation 
Note Added : Results obtained in this paper for the tail's equation of state, Wh, differ 



20 



from those reported in Lemoine et al. [21 



The fact that we include the curvature term a" /a under the definition of the generalized 
frequency u 2 , while they keep the two contributions separate, is not the only source of 



discrepancy. There are fundamental differences between the two studies. Authors of [[21 
ignore the crucial effects of the out-of-equilibrium dynamics and the breaking of Lorentz 
invariance, by the nonlinear short distance modes, when carrying out their calculation for 
Ph and pn- In their approach these effects would be contained in the dynamics of the 
vector field, described by a lagrangian C u for this field. Obviously, the choice of u M 
field dynamics and its Lagrangian C u strongly depend on the details of the short distance 
nonlinear model considered. The expression for C u these authors borrow from the Corley- 
Jacobson (CJ) model is consistent only with the dispersion function of the CJ model 
for stationary backgrounds, since both C cor and the terms given by C u contain up to 4-th 
order derivatives, together with the antisymmetric tensor F^. Therefore, it should come 
as no surprise that C u gives zero corrections to pn, Ph when applied to the FLRW Universe 
(where clearly the antisymmetric tensor vanishes identically F^) and hence, higher order 
counterterms in C u are ignored, while at the same time in C cor they have a series 

of all higher order derivative terms, up to n — > oo. It is easy to check they break energy 
conservation law and Bianchi identity, by plugging in the energy conservation equation their 
expression for p and p in their Eqs. (27) and (28) 

(p) + 3H((p + p))^0 (59) 



which we suspect is due to the abovementioned reasons. 

Perhaps they could reconstruct their formalism, either by identifying the correct C u 
which would be appropriate for the model given by their £ cor , or by defining an inner 
product (following the construction in |L5[], (0i n in ) = / dx(pF(k)4>) for the fields in such a 
way that it accounts for the nonlinear dynamics and Lorentz noninvariance of short distance 
modes with frequency F(k), in an expanding Universe. However these issues extend beyond 
the scope of this paper and we do not intend to elaborate further. 
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5 Appendix 

5.1 Distribution function for the Linear, Crossover and Trans- 
planckian regime. 

We have a time translation Killing vector for future infinity that determines our outgoing 
positive frequency modes. Let us consider our Universe as an expanding box, with size 
L = a/M, filled with modes. Inside this box we have a smaller box with fixed size l p ~ 1/M 
that determines the range for the transplanckian (TP) modes. There is a preferred frame 
attached to the small box (due to the breaking of Lorentz invariance by the short- distance 
modes) but Lorentz invariance is restored in the big box, the expanding Universe. Following 
along the arguments and derivation in [15"] for the 'particle' distribution function, it is 



straightforward to apply their expression to our model. Below we consider three regimes 
depending on the value of the momentum p with respect to the cuttoff scale M 

(a) p/M <1, the "normal regime". In this regime the frequency of the modes becomes 
linear, 

iu[p]~pe- p / M ^ P/M «i P- (60) 
The wavelength of these modes is then A ~ O(L) ^> lp. 

(b) p/M ~ 0(1), the crossover or intermediate regime during which the TP modes go 
from the "TP box" with fixed size lp = 1/Mp into the "big box" with size L = a/M. This 
process occurs due to the redshifting of the modes, Pi = ki/a. Each mode Pi will crossover 
and become a "normal mode" at some time = hi/M. 

(c) p/M > 1, the TP regime, such that A ~ 1/k < l P = 1/M. 

We do not report the derivation of the distribution function since the reader can find it 
in detail in [O . In what follows we apply it to our case, to lend support to our assumption 
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that the TP modes, modes in the small box with fixed size l p = 1/M with respect to the 
preferred frame, are not in thermal equilibrium, while modes in the range of the linear 
regime, (modes in the big box with size L=a/M, where Lorentz invariance is restored, the 
FRW Universe) are in a nearly thermal equilibrium situation. 



N{<J) 



)Vg{pv 



(Uout)Vg{Pout) 



OL, 



(61) 



The index in (out) refers to the incoming (outcoming) modes as defined in 0; v g is the 
group velocity, 

/ x did 

= (62) 

and f3 p , a p are the Bogoliubov coefficients, which in our model do not depend on the 
momentum p 0. The dispersed frequency is given in Eq. (|5]), and thus 
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'out ) — 1 ) 



since ui(p out ) — p ou t- The phase velocity is defined as v c = u/p. Therefore 

2 



N(oo) 



1% 



a. 



f e -p/Mp\i e -p/M P 



P 



The thermal distribution is immediately recovered in the limit p/M 1, i.e., v g - 
u(p) -> p. 

Now we can evaluate the distribution function for the different regimes above: 
(a) In this case 



(63) 
(64) 

(65) 
1 and 



Pp 



a, 



(66) 



as it should be since we recover in the out region the normal linear frequency for the modes, 
p/M <C 1 and v g ~ 1. 



(b) In the intermediate crossover regime, p/M ~ 0(1), we have 

A4H r 



(67) 



2 eP" - 1 ' 

where we have identified the inverse of temperature (3 ~ a and the linear term in p/M with 



M 



1 - p/M aM - k 



(5 ~ a. 



(68) 



Clearly in this regime \v g \ goes to zero, and (3u goes to infinity. The spectrum is nearly 
thermal, however M(uj) goes to zero, since as seen from the TP box the group velocity of 
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these modes as they approach p ~ M is becoming zero; or as seen from the normal particles 
in the big box, these modes have a high frequency (u ~ M) thus they do not contribute 
very much to the energy of modes in (a). 

(c) However in the TP regime the distribution function strongly deviates from that of 
thermal equilibrium, since in this case v g (in) 7^ 1, and the frequency is highly nonlinear: 



AC 



UJ 



C 
2 



\ e -p/M\ e -p/M 






M 



(69) 



with p/M ^> 1. Nevertheless M c {uj) — > when p/M — > oo, thus their contribution to the 
energy is suppressed. The suppression comes directly from the frequency u>{p) ~ pe 
this case. 



-p/M 



111 



The volume element in momentum space, dVp, for the dispersed 'particles', whose world- 
line intersect an hypersurface element rfE around x, having momenta in the range (p,p + dp) 
is dV p = 25(p fl p fl )dp 4: , where p is future directed. Based on the definition of [15] for the inner 
product of the fields 



\out"\ 
Lin 



\out 
Lin 



]) and integrating over the entire mass shell, the 3- volume 



in momentum space is given by: 



d 3 p 



(70) 



which is consistent with the quantum field theory expressions of currents and 'particle' 
number densities, given in terms of creation and annihilation operators. 

The distribution function, Eq. (|5"9"P, justifies our assumptions of deviation from thermal 
equilibrium in the TP regime. 



5.2 Averaging of the TP energies, pressure and equation of state. 

We have shown that modes in the TP regime of such very short wavelength \?p <C lp are 
out of thermal equilibrium. Thus we need to average their effective contribution to the 
energy and pressure over many wavelengths, in order to obtain a nearly thermal, large scale 
state. This averaging is done over many wavelengths since clearly scales of cosmological 
order that are of interest to us, are much much longer than any TP wavelengths. In what 
follows, we are interested to find the averaged bare quantities (p^ and (pi) before including 
any modification IT. The viscous pressure modifications,!^, occur due to the freeze-out and 
the change in the number of particles and are accounted for separately in Sect. 2. Therefore 
region II will be grouped together with region H, since both have a highly dispersed TP 
frequency and, for the moment, we are ignoring the freeze-out corrections. Approximately 
we can write the energy of regions (0+1) and (II+H) in one compact form to avoid repetition: 

p = p A Q(M -p) + M- p z e -(p/m P )e( P -M) ) (71) 

where Q(p — M) is the unit step function that takes the value Q(p — M) = 1 for p > M and 
zero otherwise. Clearly, for modes with p/M > 1 we get pn (ph)- And for modes M > p 
we get the radiation energy density of the (nearly) linear modes. 
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Using the energy conservation law (while ignoring the freeze-out effects): 

+ _ pdp _ p I dpi | dp 2 \ 
P ^ 3 dp 3 I dp dp J 



with 



= (Pi + Pi) + (P2 + P2) , (72) 

Pi = fpV^, p>M, (73) 
P2 = p\ p<M. (74) 



And from Eq. (|72|) we find: 



wi = ^ = -^Q(k B -aM), (75) 
pi 3aM 

w 2 = — = -. (76) 
P2 3 

It is clear from the previous section that since the (nearly) linear modes, p 2 are nearly in 
thermal equilibrium then we do not need to bother with the averaging procedure for them. 
(One can however check to verify the result (w 2 ) = 1/3). This is not the case for the modes 
with p/M > 1, since these short wavelengths are out of equilibrium. Let us calculate (pi), 
(pi), and (wi): 

, _ g Pl a 2 da _ J a * =k/M Pl a 2 da 14M 4 q* 3 

\Pl) - ra oi„ - 6 ~i - 1C „„Q ' {<<) 



Jq a 2 da a 3 16ea 3 



Similarly, 



and 



/_ v IoPa 2 da Jq h/M pa 2 da 14a* 

Pi = ra 2 , =3 § =— ^-(M- 7 § 

J a z da a 6 20a 

W = ^ = 0-^.-5, (79) 
(pi) 2Ua a 

where in Planck units, a* = l p and we have used the normalization that Planck length 
l p = 1/M = a* = a(tp) = 1. The time a* corresponds to the crossover time when the 
wavelength of the TP modes is on average of order the size of the small box, p(a*) = M. 
This time scale of order the Planck box is much smaller than scales of cosmo logical interest, 
L = a/M. External observers bound to the large scale, Lorentz invariant Universe with 
size L = a/M and in thermal equilibrium, 'feel' the energy and pressure contributions from 
the TP modes given by Eq. (|69|). 

Of course the real equation of state for these modes is given by their effective equation 
of state, Wi, Sect. 2, that takes into account the relativistic kinetic theory modifications due 
to the freeze-out. 

In a similar manner, one can obtain (w 2 ) and in particular the numerical coefficient A 
for (w/),p/,p/. 
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